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USE 01!' FROHLICH-RAIMES AND OTHER ELECTRONIC MODELS 
FOR PHYSICAL PROPERTIES OF METALS AND ALLOYS''.' 


James T. Waber 


+ 


The offecliiveness of the Frohlicb-Raimes model in treating 
force constants and alloy behavior is discussed. Typical predic- 
tions of trends which have been made are; of the deviation of 
Vegard's Law and the variation of the bulk modulus with pressure 
and with alloying elements. The simple model of Koskimaki and 
Waber using the linear combination of the density of states has 
lead to several useful predictions for titanium-based alloys. 

This model will be discussed and recent results are reviewed. 


* This research was supported in part by the National Aeronautics 
and Space Administration and by the U. S, Atomic Energy 
Conmission. The work was initiated while the author was at 
Los Alamos Scientific Laboratory, 
f Department of Materials Science, Northwestern University, 
Evanston, II. 60201 
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Several years ago Brooks and Wigner and ScitE took up an 
investigation on the effectiveness of a volume -dependent theory 

3 

of metals originally developed by Frohlich in 1937 and modified 
4 5 

by Raimes ’ to include several conduction electrons per atom. 

6 

In an unpublished thesis, Bernstein applied this Frohlich- 
Raimes (FR) method to rare earth metals. Brooks notes that the 
equilibrium lattice spacing in normal metals is essentially 
determined by pressure of the Fermi gas of free electrons. 

•7 

In 1961, Larson aiid Waber reinvestigated the problem and 
noted ,ui undetected asyiriptotic singularity existed in the basic 
equation. Wlien this was removed by retaining more terms 
before truncation of a series expansion for the wave function, 
better agreement with experimental data was obtained. 

8 

Another modification by Waber, Larson and Smith was to der- 
ive the equations assuming that the state at tlie bottom of the 
conduction band could have an angular momentum greater than 1=0, 
This permitted them to investigate transition metals. A brief 
review of the theory and some results will be given below before 
applying it to alloys, since the recent developments have not 
been reported in the literature. 

Before continuing with that topic, it is desirable to give 
a little background on the second model which will be used inas- 
much as it is apparently unknown except to a small circle of 
colleagues. In the case of metals, the Pauli exclusion prin- 
ciple requires that a distinct sot of quantum numljcrs be assigned 
for each conduction electron coming from an atom in the solid. 
Instead of the usual set of n. A, m and spin for free atoms, they 
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arc k , It , k and spin, and arc associated with the momentum k 

x’ y’ 55 ’ 

of itinerant electrons. TC the density o£ atoms per cubic cen- 
timeter oecomes larp.c enough* the discrete spectrum of energies asso- 
cia.i;ed with each k value can be replaced by, a piecewise continuous 
E(k) curve, Tlie number of discrete states within any energy 
range E to E+dE is replaced by the density of states function 


N(E). The band structure of metals is concerned with calculating 
such N(E) curves for a pure metal and with deducing various 
observable quantities from them. 


When the original proponents of this theory, such as Jones 
10 

,1 ai^d Hume-l’othery theory applied it to metallic alloys, they assumed 

11 

a Rigid Band model (RB) . Sinter made onp of the first uses of 

this model to explain the variation of magnetism of alloys; the 

12 

famous Slater-Pauling curve wos the result. These men were 
keenly aware of the limitations of the model, but insufficient 
information was available for them to regard it initially as any- 
thing but a working hypothesis. It has proved to be very success- 
ful in making important predictions during the last thirty-five 
years and many experimentalists have verified the essential 
features; the comments of Raimes are interesting. However, 
repeatedly small and sometimes large deviations from the Rigid 
Band prediction have been noted. 


In an extensive series of band calculations for 3d transition 
metals, Waber and Snow'^ ^ showed that the band structure was far 
from rigid for such a series of metals; but they also showed 
why the Rigid Band (RB) model worked well. Despite changes in 
the shape and the location of the bands for individual metals, the 
plot of N(Ep) at the Fermi level E^, when plotted against the 
atomic number had the two hump shape of the origiiial E(k) curve 
computed for pure metals by Manning^^ and by Chodorow^^. 


Waber proposed that it would be preferable to use a linear 
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combination o£ the two (or three) density of states for the 
constituents in dealing with alloys. In discussions during the 
conference on the Electronic Density of States held at the 

17 

Bureau of Standards, Waber used the LCDS model to estimate (a) den- 
sities of states curves which have three humps and (b) a typical 
curve of N(Ep) versus composition for one type of titanium alloy. I 

He also presented the three peak N(E) curves for several 4d metals, 

18 

Koskimaki and Waber have computed these for a number of alloy j 

19 i 

systems. Colli, ags and Gegel have made effective use of the j 

calculations recently. It is appropriate to present further the j 

justification of this LCDS model and to record some of the results | 

obtained with it. Independently, Sterns^^ studied a form of the LCDS | 

model and concluded that should bo constant. | 


Let us contrast these two approaches to the behavior of alloys. 
They depend very directly on the average spatial density of elec- 
trons in the metal. The modifications of the Frohlich-Raimes 

Theory ignore all ciuestions about the arrangement of atoms in 

2 2 

the solid. It is also assumed that the Fermi energy (h k p) 
depends only on the number of electrons per atom (i,e., on the 
(c/a) ratio). 




! 


* 


The LCDS theory is predicated on the detailed E(k) curves 

for each constituent, which in turn is based on the local ar- 

14 

rangement of atoms. That is. Snow and Waber have shown that 
the E(k) curve is not the same for the BCC and FCC forms of the 

I 

same metal. S 

f 

Thus, the LCDS model is better able to cope with alloys ani f 

' I 

such phenomena as short range order and clustering, than is the P 

■volume dependent FR model. It in turn is more effective than the I 

(e/a) dependent RB model. Each model has a range of problems | 

for which it is effective. For the calculation of cohesive j= 

energy, compressibility, etc., the FR approximation will be 
very useful. For the behavior of the specific heat and magnetism 
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o£ alloys, the LCDS model will be a distinct improvement over the 

RB model. In principle, the volume (or pressure) dcpendcuW oC li(k) 

curves can bo calculated by a band calculation such as by the APW 
21 

method. Averill and independently Perrot have calculated 

the cohesive energy for alkali metals with considerable success. 

23 

Very recently, Snow has done it for copper. However, these were 
extensive calculations which consumed large amounts of time on 
the computer. It is not practical to do exploratory research on 
alloys by such methods. 

The two approximations outlined in contrast are very simple 
to carry out and therefore lend themselves to treating compli- 
cated problems. It is with these justifications and limitations 
in mind that we undertake to study the behavior of alloys. 

Formul atlon of Frohlieh-Raimes Model 


The basic assumption is that the behavior of the conduction 
electrons in a metal closely approximates those confined to a 
sphere of radius R„. Next one assumes that the ion core with 
its charge of No (or N, in the atomic units where e4i=m=l) does 
not occupy a significant fraction of the sphere. One solves the 
Schroedinger equation for one electron confined to the Wigner- 
Seitz sphere, which also contains an ion core. We will call 
'^.eo (\) the radial part of the wave function in this field 
(averaged over the 4 tt solid angle). Tlie subscript zero indicated 
the ground state for the angular moment f . The Schroedinger 
equation can be written as 


2 ^^fo 

,2 r dr 
dr 


r 


f, = 0 
io 


( 1 ) 
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The usual Wigner-Seitz boun<loi 7 couditlon^^*^ is applied, namely 
that the gradient vanishes 



B 0 

j’^“\ 


( 2 ) 


In the space available, it is not practical to discuss the 
mathematical part of this problem. Only the results in the form 
of the principal equations will bo indicated. 


25 

As Wabor, Larson and Smith show, the lowest eigen-value 
depends on in the following simple way 


E (R ) 
0^ s"^ 







(3) 


The parameter R^ corresponds the radius of which the energy 
becomes a minimum. At the other oxf^ome end, when R becomes 
infinitely large, approaches I^ which is the n'^ ionization 
potential for removing all but one of the valence electrons from 
the atom in question. 


2 

Let us call p the positive energy difference between the 
free ion and the E for the metallic case, and call the ioni- 
zation potential w . After substitutiiig R^ for R^ in Equation 
3, we find that 





(4) 


m i a-n:> 2 

R ~ 2 ~ 

0 


(Aa) 


6 


Waber 




4 


In the procedure outlined initially by Prohlieh, one expresses 
the eigeia-function where 

J 

P V C r“^' 

“ exp (-wr) r^ L k 

ic«0 (5) 

in terms o£ a power series P » n/w. 

The recursion formula for the coefficients in the power series 
can be shown to be 


■2w(k+l) 


Tlren one requires that the function be well behaved at the 
special value of the sphere boundary. Thus, as a generalization 
of the solution which Frohlich and Raimes found, one can write 


R (B-A) cot bR -(H-pR A) 
o o o 

R (1+Ar cot ) 
o' o o' 


dr r=R 


where A is an arbitrary small constant (which*’ it will be convenient 
, 25 

to set to zero). 


Substitution of the wave function alternatively leads to 

. J 


(1 to 




numerically by iteration, 
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By Gxpurienec one xioed ineludiv au most six tonss and in many cases 
two will suffice. These equations reduce to the ones given by Ralines 
when one sots £-0. 


A plot of the ground state radius is presented in Fig. 1, 
as a function of A discontinuity is found because the cotangent 

function approaches (i®) when its argument equals rr. This behavior 
of the logarithmic derivative corresponds to f^^ vanishing at the 
cell boundary for a specific energy. If one desires, this difficulty 
can be eliminated by retaining the terms multiplied by A in the 
logarithmic derivative, as the smc,''th curves show. Insets for other 
values of N are shown, 


The next step in this volume dependent theory is to put n~l free 
electrons back into the spherical cell. Thur'' are several types of 
interactions bctweeii the electrons. The first is Coulomb energy of 
the N electrons. The second is the mean Fermi energy of these elec- 
trons confined to sphere of radius R . The third is the free electron 
exchange energy. A small correction is also included for correlation. 
Thus, the total energy for the N valence electrons is given by the 
simple formula 


E (R ) 

b 


“ E 


(R ) 
' s 


1.2N , 2.21 


11 


+ 


r . 


0.916 

r.. 


corr 


( 9 ) 


Inherent in this expression is the a.dea that is approximately 

constant throughout a significant fraction of the sphere (for cases 
Raimes^ studied, the fraction was approximately 90 percent.) In 
uniform electron gas, we can set r^ as the radius of the Fermi hole, 
namely R^/^Fn. The simplest form of correlation was used, namely 
that based onWigner's approximation formula 


corr 


■0.288 

R -hS.l 
s 


( 10 ) 
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Willie improved trcatmenUo of exchange are available si.\ch aa 

26 27 

choGo by Slaccr and hie oolleaguea and by biborman , 

Gigniflcant precision should not. be expeebed with this over- 
simplified model. Typical interaction terms which are independent 
of the element involved, are presented in i-^ig, 2 for N«4, 

Wliilc ii: might be desirable to employ an lmp:f;/vcd 

, 28 

treatment of correlation done by Pines and 'slwres , it is 
not likely to lead to significantly different results. Just 
as Slater has indicated in his recent studies that the exchange 
term should bo multiplied by a constant one could 
correct the kinetic energy by dividing by the effective mass, 
m ^ each of the coefficients might be regarded as adjustable 
parameters. 

It has been decidad to apply the more stringent test to 
the model, namely, the original free-electron coefficients have 
been used in this heuristic study. 

The energies E^Ckg) and E^.(kg) arc presented in Fig. I for 
titanium and rirconium. The equilibrium radius corresponding to 
a minimum in is designated p. Even though the inter- 

action energies depend only on N and arc independent of the 
specific metal being studied, the curves of E^(R^) may be 
different for two very similar metals as titanium and zirconium. 

Estimation of Other Physical Quantities 

Brooks^' noted that the experimental cohesive energy per valence 
electron is surprisingly constaxit. For 23 elements it ranges 
from 0.060 for Cesium to 0.137 Rydbergs per atom for carbon, 

llie cohesive energy data tabulated for the lanthanides by 

29 30 

Gsehneidner and the data of Trulson, Hudson and Spedding 

for the heat of sublimation at 298® K were divided by the number 
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DeUL-rmination of H ^ and p from a 
of titanium and zircon 








o£ valence electrons. This energy was found to vary for the 
lanthanides from 0.106 Rydberg for lanthanum and 0.0636 for 
thulium. 

The cohesive energy S, i.e. the heat of sublimation of O^’K, 
can be obtained from the total energy at p by means of the 
exprof^sion, 

n 

S = I I. - NE,p(p) (11) 

j 

wherein the several ionization potentials I. are summed. 

•i 

The pres sure -volume relation which leads to the compres- 
sibility X at any pressure can be obtained in a straightforward 

manner from the E,„(R ) curve, Raimes^*' points out that the 
X s 

pressure P corresponding to a relative change in volume can be 
obtained from the expressions 



Note that 1 atomic unit of pressure used in (12) is equal 

to 150 X 10^ kg/cm^. No attempt has been made here to correct 

these data to absolute zero although generally speaking this 

change would tend to improve the agreement between calculated 

and experimental values. The compressibility x zero pressure 

can be obtained from the second derivative of E^(Rg). With 

one exception, the listed values of the experimental corapreesib- 

31 

ility were all taken from Bridgmen's tabulation of the 

quantity x using the more recent data for 75°C. The value for 

32 

silicon is from Pearson and Brattain 

Experimental and calculated values of the three physical 

10 
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coivoCants p, S and x prescnund ii\ Table 1 lor a number of 
metals and si'ml-metals. In as much as S is the di££erovaoa o£ 
two large terms, the cflect ol any error In estimating the indi- 
vidual electronic energies is greatly magniiied. To partially 
overcome tills objection, an experimental value ol 12,., was obtained 
by adding the several experimental ionisation potentials to 
the experimental heat ol sublimation. These two values of 
E,j,(p) are included in the table. The energy conversion factor 
used was one Rydberg per atom equals 313.66 kcal/mole. The 

values of the heat of sublimation S for monatomic gaseous spe- 

33 

cies and co'rrectcd to 298® K were taken from Teatum et al . 


.iWliercver possible, the ionixation potentials obtained by C. E. Hoore- 

Sittorly^^^ from spectroscopic data were used. One of the most 

readily usable and complete tabulations of ionization potentials 

35 

lias been presented by Einkclnburg and coworkers . Tlic ioniz- 


ation pr.tentials for Hf m and Ilf IV were recently given by 

36 

Klinkcnberg, van Kleef and Noormnn . The ionization potentials 
for thorinm are not well known, the sum is probably accurate to 
bettor than 5 per cent. 


Values of V, S and p were obtained independently by Raimes 

A 

who used the simplified theory cited above . For comparison, these 
values for trivalent and tetravalent metals are presented in Table H. 


In the calculation by Bernstein^, of the cohesive energy of 
the lanthanides, the experimental values of the third ionization 
potential I„ were used, where available^ to obtain the R . Waber and 
Larson recalculated p, x ^tid R^ for a number of rare earth metals. 
These data are presented in Table III . Values are given for both di- 
valent and trivalent forms of ytterbium and europium. The use of both 
valences is much more questionable in the case of luteteum. 


5 

Raimes pointed out that, aside from the systematic errors 
due to the approximations made in deriving the equations, one 
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'i'ublc 1. Comparison vU’ the Ualaulateil Valuos of lUe 

Motallle Railiur. p, lloai of b'ubiimuuion B 

ami CompresuiibUUy \ wi th the Kxperiwcntuil 

Valutiti. Valutn) of ibu *rotal lilac trotiic 

Energy K,, ami tlrouml Biauc Rmllua r also 
* 0 

UiK'tl, 


Klemenl 

Radltui Radius* HeaU B 

Oompreaii,** 

ToUal R.j, 
kCiil/raole 

vJousl tiered 

A,U, p A.U. kCal/mole 

\ X lO^’cm^VUg* 


Trlva 1 iMU~ Kl ement a 


Aluminum 
Exp , 

?.13R 

3,1')Q 

2.900 

26 

75 

0.78 

1.34 

1253 

1302 

Beandlum 

Kxp. 

. bhO 

3.515 

3.427 

ino 
80 . a 

1.27 

111? 

109? 

(la Ilium 
Hxp . 

I.DJ3 

<Mt •* M 

2.973 

3.1U1 

13 

66 

0. 61 
2.0 

1332 

138'4 

Yciritim 
Kxp . 

3.1bB 

«■ M M «« «* 

4.358 
3 . 760 

47.7 

82,6 

2.44. 

2.95 

949.2 

984,1 

Lt\vlium 
Kxp , 

3.173 

*«•••»•««» 

3.215 

3.47H 

24 

59,1 

0.82 
2, 5 

1238 

1273 

Lauihanum 

Kxp, 

3,466 

4.653 

3.9.11 

60 

99.5 

3.22 

3.24 

894 

922 

'rival Hum 
Exp , 

2. DOB 

3.056 

3.577 

3 

43 

0.66 

2.77 

1301 

1341 




Tetravalenf HI emeu la 


Silicon 

Exp. 

1.846 

mi. am 

2.875 -46 0.298 

3.184 105 1.00 

2354.8 

2497.0 
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Table I. (eonl.inuo.tl) 


Element 

Cousidored 

Radius 

r 1 A , d • 
o’ 

UaUivu? 

p A • • 

lie at: 8 
kCai/molo 

Compress, 

X X 10 cm /kg 

Total E,j, 
kCal/mole 

Tiunlum 

l.‘)4d 

2.98B 

95 

0.342 

2195 

Exp . 


3.178 

Ul.l 

0,868 

OQ1 •) 

•a M iM 

Cleirmaniim 

1.B16 

2.842 

-21 

0.290 

2370 

Exp . 

----- 

3.308 

89 

1.364 

2480 

Zirconium 

3. bis 

3.735 

82 

’0.810 

1864 

Exp . 

-- — 

3.347 

142.2 

1 . 106 

1942 

lla Cnium 

.3,072 

3.799 

29 

0.816 

1836 

Exp . 



3.301 

168 

0.881 

1975 

Thorium 

3,077 

4.2.56 

17 

1.321 

1657 

Exp. 


3. 754 

127 

1.846 

(1760) 



Pontavalont Klomonts * 


Vanadium 

1.634 

2.628 

169 

0.14 

3869 

Exp, 


2.817 

122 

0.609 

3892 

Niobium 

2.237 

3.310 

75 

0.34 

3163 

Exp . 


3.071 

185 

0.57 

3273 

‘V Tha expo 

rimo,nUal values o£ 

the Wignor- 

SeiLu radius were 

obtained fro 

the volumes V cm 

per gram 

atom lii'.Lod 

by TeaLum, Otjchtieidncr and 

Wabar , 

using the 

formula 

p«^1.3880 



>'“'VExporimonl:al compi, os s ib i 1 i ties 


are ior 343'’ K, 
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Table II Unpublished Calculations by Uaimcs 


METAbS OF IIICHKR VALKNCY 


Me till 

Lattice Constant 

p(S> 

Compressibility 
10^“X (cm^/dyi\e) 

Cohesive Energy 
S(Kcal./molc) 

Deviation 
AS as % 
Fermi E 


Calc. 

Ob fi . 

Calc. 

Ohs. 

Calc. 

Obs . 


Trivolent Metals 







A1 

1.72 

1.58 

0.9 

1.4 

8 

55 

-12 

Sc 

1.99 

(1.8) 

1.5 

— 

75 

70 

2 

Y 

2.36 

1.97 

2.8 

— 

28 

90 

-29 

Ga 

1.60 

1.67 

0.7 

2.1 

- 6 

52 

-12 

In 

1.74 

1.84 

0.9 

2.7 

6 

52 

-11 

T1 

1.65 

1.89 

0.7 

2.3 

- 17 

40 

-13 

Tetrnvalent Metals 






Ti 

1 .61 

1.62 

0.4 

— 

132 

100 

4 

Zt 

2.02 

1.71 

0.9 

— 

14 

110 

-20 

Sn 

1.74 

1.86 

0,5 

1.9 

- 36 

78 

-18 

V 

1.42 

1.48 

0.2 

— 

100 

H 

85 

+ 1 

Semi- 

•Conductors 







Si 

1.55 

1.68 

0.3 

0.3 

- 63 

85 

-18 

Go 

1.53 

1.74 

0.3 

1.4 

- 59 

85 

-17 


*The author is indebted to TroC. Raimes who made this information 
available for use. 
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TabXcin. Comparison of Results Obtained in Recent Cohesive 
Energy Calculations for the Lanthanide Metals 


Wigner-Seitz Radius p Compressibility Parameter R^ 



(Bohr Units) 



10^ cm 


(Bohr Units) 


BS 

LW 

OBS. 

BS 

LW 

OBS. 

LW 

BS 

Sc(3) 

3.69 

3.635 

3.427 

1.37 

1.27 

— 

2.54 

2.56 

Y(3) 

4.40 

4.358 

3.760 

2.69 

2.44 

2.95 

3.198 

3.20 

La(3) 

4.70 

4.653 

3.921 

3.47 

3.22 

4.13 

3.466 

3.47 

Ce(3) 

4.65 

4.464 

3.81(y) 

3.29 

2.75’ 

4. 10 (y) 3.295 

3.42 

Ce(4) 

4.31 

3.490 

3.57 (cv) 

1.58 

0.61’ 

4.72(cv) 2.387 

3,09 

Sm(3) 

4.50 

4.174 

3.761 

2.88 

2.15 

3.52 

3.684 

3,28 

Eu(2) 

4.84 

4.829 

4.290 

8.13 

7.76 

4.91 

3.715 

3.67 

Eu(3) 


4.088 

4.290 

— 

1.98 

4.91 

2.953 

— 

Yb(2) 

4.54 

4.485 

4.044 

6.33 

5.92 

7.84 

3.389 

3.28 

Vb(3) 

— 

3.986 

3.986 

— 

1.80 

7.84 

2.860 

— 

Lu(3) 

4.25 

3.771 

3.628 

2.36 

2.51 

2,33 

3.198 

3.06 

Lu(2) 

— 

4.358 

3,628 

— 

3.10 

2.33 

2.701 



BS K reference 6 
LW = reference 38 
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soul'co ol’ error in the coheaive energy ariucs from the eritimutlon 
of the i, round au.iLe enerfiy; a two per cent ehaniio in K^Cp) 

■ yields almost a 20 per cent change in S. It is shown in Fig. 1 

that for certain regions of the ioniaation potential, an appro- 
etablc e.rror could be made in determining by means of the 
Raimos equation due to the steps. Such errors would have a 
stronger effect on the ground state energy at R “ R , than at 
R » p . However, hero wo will not go into problems of estimating 
R^ accurately. 

Brooks^'"^ presented a simplified theory of cohesion which 

, involved a modification of the quantum defect method used by 

Kuhi/^^ to determine E (R ) . Values of r determined in this 

0 0 o » 

m\y are compared in Table IV with those computed herein. The 
‘i equilibrium radius p was evaluated using Equations (8) and (9) 

as above. The values of p obtained in these two ways are also 
’ compared with the experimental value in Table IV. It will bo seen 

that Brooks' small r^ values result in bettor agreement with the 
experimental radii. Also Lite experimental values of S for the 

* five metals studied by Brooks agree with his calculated values 

.1 

^ better than they do with those calculated by the present method. 

, The Wigncr-Scita radius p and cohesive energy of the noble 

metals, copper, silver aiid gold are compared with experimental 

1 

j values in T able V. Brooks noted that the experimental values 

’ of S for these raonovulont metals were three or four times larger 

• than the values of S/N for any of the other 32 elements he con- 

j 

5 sidered. lie concluded that the closed d-shells contribute sig- 

i nificantly to the binding. It is interesting that for the 

i IIB, IIIB and IVB elements which follow these transition mofals 

I . in the Beriodic Table, the binding is normal and the values of 

I S/N lie in the same range as mentioned in the first part of Lliia 

1 , paper. Thus the ion cores are more effectively separated from 

# 

1 each othaj: by t:ha gas of free electrons; an error arises when 


j 


\ 

I 

I 

i 


! 
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Table IV, Comparison oJ; the brooks' Values o£ and p with 

Present Computations of these Radii in Atomic Units 


Element 

Treated 

Ground Sl;»te Radius R 

0 

Brooks Preseni 

Wip.ner-Soitz Radius, p 
Present Brooks Exper, 

Na 

2.97 

3.244 

4.147 

3.886 

3.991 

Mb 

2.38 

2.583 

3.601 

3.437 

3,345 

A1 

2.03 

2.132 

3.190 

3.088 

2.990 

Si 

1.755 

1.846 

2.875 

2.75 

3.184 

Zn 


2.003 

3.024 

2.784 

2.904 
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Table V. Comparioon oi' CaleulaCed and Experimental 

1 

Valuoa Tor the Noble Motala (after Brooks ) 


Cohesive Energy, S Wlp,nor«Scit:; Compressibility 


RydberBs/olecuron 

Radius, Bolir Units 

X X 10 cin /kp, 

Calc, Exp. 

Calc. Exp. 

Calc, Exp. 


Copper 

0.087 

0.260 

2.608 

2.669 

2.80 

0.734 

Silver 

0.086 

0.221 

2.666 

. 3.017 

3.03 

1.004 

Cold 

0.085 

0.264 

2.146 

3.012 

1.43 

0.570 
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one asiijuiuojj that; the* <,l“cXactrortf> avo conCincJ to a vci*y ywull vcfiioft 
assoeiateU with th« ion core and have nci’.lihihlo probability o£ being 
£ound at the cell boundary. This io equivalent to the old tight 
binding approximation which ia inadequate tor transition metals. 

Typical behavior of the calculated compressibility \ as a tunction 
o£ pressure is shown in Fig. 4 Cor Cour tetravalcnt metals. 


in the noble metals, the interatomic distance, however, is gov- 
erned by the interaction between d -electrons and the ion cores. Con- 
trary to the treiid found Cor most metals in Table I, the calculated 
values of p in Tuble V are less than the observed values. Thus, one 
would anticipate that, duo to ionic repulsion, the actual compressibil- 
ity of copper, silver and gold would be su})Stantinlly lower than that 
calculated by the present method, which only includes electron-electron 
repulsicns. This is borne out in Table V where the cxperlmontal x 
values are only 0.3 to 0.4 that of the calculated values. Thus, it is 
roasonatle to conclude that for these noble metals, the interaction 
between the valence electrons and those forming the ionic core is much 
stronger than in normal metals and that the present .“implo method is 
not satisfactory. 

The importance of such repulsion is illustrated in the pressure- 
volume graph, Fig. 3, for nickel. As a shock wave passes through a 
metal, it is subjected to large negative pressures during the rarefac- 
tion wave --although that portioii is itot plotted in Fig. 5. As ai\ 
atom, nickel has only one 4s electron to act as a valence electron. To 
incorporate the effect of the 3d electrons, wo have assumed that N is 

higher than two. In doing the calculations, N®3 gives better agreement 

41 

with the experimental shock wave data . 


A more significant way of dealing with the ion cores is to add 


K. . 

XOlliC 


b exp 

to the terms in Equatioii 9. 



( 14 ) 
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V?aher 


Variation of C ompressibiHt 



Pressure 






ro 

o 6 


Here b io an adjua table eonntanfc approximately equal to 10^^ ergo 
and S.| . is given by 


where r and r, ore the ionic radii euch ao moot chemists use, Con- 
■i* j 

coming the repulalve parameter which occurs In the Bom-Mayer 
expression for ionic crystals, llafnoister and Zahrt*^^ showed that 
values of pjjj,j could be calculated directly from overlap integrals -- 
nomtally It is an experimentally adjustable quantity approximately 
equal to 0.345A®. Wc have established that Including ionic terms im- 
proves the agreement with shock wave data on gold as Fig. 6 shows. 
Inasmuch as the inclusion of ionic repulsion has been made only for 
limited cases here, it would take us away^from the main theme of this 
paper to go into further detail about how the quantities b and 
are estimated, I'his will bo discussed in a forthcoming paper. 


/.p olication of the Ffdhlich-Raimes Model to Alloys 


38 

The Frohlich-Raimes method has been applied to alloys . One 
reasonable assumption is that the ground state energy of the alloy is 
given by a linear combination of terms, Th*at is, 


-E„ (CM 


C n. 
A A 

II 


R‘ 


k 


oA 



(16) 


where C^, n^ and art, respectively, the atomic concentration, the 
valence and the Frb'hlich-Ralmos parameter for component A, and is 
the variable radius of a Wigner-Soitz sphere. Here it is assumed that 
i=0. For the interaction energy, the mean number C^n^ + of 
valence olectrons wa.s used. This n is inserted in the remaining terms 
of Equation 9 rather than averaging the two E^, values. To distinguish 
expressions, we will denote the total energy of an alloy by Values 
E , radius p, and compressibility y have, been obtained as a function 
of concentration for a large number of alloys. Most of these results 
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will bo reported elsewhere, but one type of result will be suinmartzcd 
here. 

It has become common for metallurgists to compare the observed 
average lattice spacing of alloys with the spacing of alloys with the 
spacing obtained by averaging the lattice spacings of the components. 

One can calculate p(C^) and compare it with Vegard's Law which states 
that the volumes of components are additive. We have assessed this 
deviation from Vegard's Law by computing 

-y(V -P(c^> - (c^p^ + (17) 

Such deviations, Ap(C^) are illustrated in Fig. 7 for Ti and two dif- 
ferent valences of cerium. A large negative deviation can be seen for 
the tctravalent Cc. 

♦ 

The total energy E^. (C^,p) has been plotted for three other binary 
systems. Pig, 8 is for Na-Ag for which little solubility exists at 
either end. However, the calculated curve cannot illustrate this fea- 
ture, The large positive deviation from Vegard's Law which can be seen 
in the plot of p (C^) is the clue. Of course, strain energy in the lat- 
tice has not been included in the present model. So direct evidence 
for phase separation cannot come from the intentionally restricted model 
used at the present time. Fig. 9 is for Na-Mg alloys. The heat of 
formation of these alloys /\E is the deviation from the straight line 
which would be given by a mechanical mixture. In almost all cases 
studied, the solid solutions would tend to be more stable than a two 
phase mixture. Note that this stability does not result from including 
the usual entropy of mixing but is strictly electronic in origin. No 
electronic entropy term has even been included. Both would increase 
the stability of any solid solution. 

The results of a large number of calculations about Vegard's Law 
are summarized in Table VI. In some coses, there is experimental 
evidence for positive and negative deviations. In a number of cases 
predictions are made which may be verified using additional lattice 
parameter data not readily available to the author when this summary 
was originally made. 
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Table VI. 


on o£ Calculated Deviation from VGRard's Law with Observed Deviation 


C. at Max. 
Calc. Ap3 


Gs-Rb 
Na Mg 
Na Ag 
Na Zn 


Mg -Li 
Mg Zn 
j Mg Sn (2) 
i Mg A1 
; Mg Ga , 

: Mg In 
■ Mg T1 

i 

i Mg La 
j Mg Ce (3) 
! Mg Ce (4) 
I Mg Sn (4) 
1 Mg Bi 


0.55 Na 
0.45 Na 
0.55 Na 


0J45 Mg 


0.55 Mg 
0.55 Mg 
0,55 Mg 
0.55 Mg 
0„5 

0,55 Mg 
0,55 Mg 
0.55 Mg 


Calc. 

100 4[p(C,)3 


+0.15 


+1.08 

+0.43 


-3.38 

+1.65 

-0.04 

-27.00 

-3.49 

-2.53 

-3.27 

+12.80 

+10.20 

-1.75 


-3.39 


Exnerirr.antal Dv*^vi 


Cubic Hexagonal 

a^ ale 

o o 0 


Monovalent Solvent 
neg. j I 


Calc.., 

100 aCp"^(cJ] 


-0.58 

-0.24 


Divalent Solvent 


Deviation r. 


No 

Vol. I Cu 


only one point 
only one point 


-1.73 

-6.98 


- 0.01 

-1.38 

- 2.11 


-1.83 
+ 3 . 3 2 
+2.77 


-2.46 

-1.97 


Ca Sr 
Ca Ba 
Ba Sr 
Ca Zn 
Ca Ti 
Ca Mn 
Ca Cr« 


0.55 Ca 
0.55 Ca 
0.55 Ca 


40.40 

+2.50 

+0.91 

+9.85 

-13.09 

-21.26 

-30.51 


-0.17 

-1.23 

-0,48 

-2.40 

-12.67 

-20.76 

-23.39 


3 3 

Note;' The linear combination of volumes C p (C.) + C p (C„) does' not lead to the same 

A A B B 

deviations (occasionally even differing in sign) as does the linear combination 
of cell diameters as used in Eq. 
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Tabic VI. (cent.) 


Comparison o£ Calculato.d Davlablon £rom Vt»>nrtl*s Law with Obsorved Dcviaslon 


1 3 in ary 
1 Sys tc\r> 

1 

V. ac 
Calc. 


[■■iSlSIlBHBBi 




Calc.„ 
100 aCp"" 


\mmsm 



IHMSIil 

pi 

Cell 

Vol. 

*SOA“* 

mC 










3 




TrivalenU Solvent: 





i 

|a 1 Ca 

0.45 

Al 

-4.09 

nog. 

— 

— 

-7.38 


■ 


A1 Zn 

0,45 

Al 

+1.59 

pos. 

— 

— 

+0.40 




Al-La 

0.50 


+8.36 

insi 

ifficicn 

, data 

-2.35 




Al-Ce (3) 

0.50 


+6.55 

pos. 

— 


-1.73 




Al-Ce (4) 

0.50 


+2.25 

pos. 


— 

+0.54 




Al-Ti 

0.55 

Al 

-1.09 

none 

neg. 

pos. 

-0.41 




Al Mn 

0J55 

Al 

-5.64 

— 

— 

— 

-2.92 




A1 Cr 

0.55 

Al 

+9.19 

— 


---- 

-3.52 









Tetraj 

?onal 










a 

c 










o 

0 





In-Cd 

0.45 

In 

-1.05 

— 

pos. 

— 

-0.39 




In-Sn (2) 

0.45 

In 

-3.01 


pos. 

pos. 

-1.65 




In-Pb (2) 

0.45 

In 

+2.54 

— 

pos. 

pos. 

-1.24 




Iii-Tl 

0.50 


+0.17 

— 

pos. 

pos. 

- -0.03 




In-Sn (4) 

0.50 


-0.49 

— 

pos. 

pos. 

-0.16 




In-Pb (4) 

0.50 


-1.02 

— 

pos. 

pos . 

-0.38 




In-Sb 

0.55 

In 

-2.52 

— 

— 

— 

-0.87 








Tetravalent 

Solvent 





Ti Ca 

0.45 

Ti 

-13.09 

— 

neg. 

neg. 

-12.67 




Ti Zn 

0.42 

Ti 

-0.53 

ins 

afficien 

t data 

-0.15 




Ti-Sn (2) 

0.40 

Ti 

— 

— 

neg. 

pos. 

-4.04 




Ti Al 

0.45 

Ti 

-1.09 

— 

neg. 

pos. 

-0.41 




Ti^-Ce (3) 

0.45 

Ti 

— 

ins 

uf ficien 

t data 

-6.22 




Ti La 

0.45 

Ti 

— 

— 

— 

— 

-7.66 




Ti Zx 

0.50 


+2.60 

none 

none 

none 

-0.52 
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n.n,n.nw.„n of Cnlclato.l ...wlncion Cr.» Vosnrd' r . La» «Kh Ol'SC-rvcJ DovlaClon 


G ”7t: >:u}j7 1 Calc; FKncriT^cnjjlJV^ ^alc.^ 

Calc, 300 ^CpCC^)] Cubic Uc;:a>oaai J-00 ^Cp (C^)J 


Dcviaciwii 



0.50 
0.50 
0.50 
0.50 
0.50 
0.55 Ti 
O.'ijO 
0.55 Ti 
0.55 Ti 


+3.03 


-0.63 
+2. 16 
+3.57 


^ I 
a 


TctravalcnC Solvent: (conU.) 
insufficient* 'data 1 


neg. nag. 


neg. nag. 


none none 



neg. nag. 


- 0.21 

- 0.22 

-0.62 

-1.64 

-0.54 

- 1.06 

-0.40 


'1-0. 1C 


2r Ce (3) 
Zv Ce (4) 
Zr Th 


0.50 


0.50 


0.45 Zr 
0.45 Zr 


+ 0.02 

+1.36 


-O.OC 

= 0.95 

-l.l! 

- 4.01 


o « 





The caleulaued data in Table Vl are ai*van«ed by the valence of the 
r.olvenU, The caaos treated were tor binary ayatems cited by I’earaon 
iw his large tabulation of lattice iiaraincters. Without becoming involved 
in detailed evaluation of the experimental deviations of different al- 
loys, as well as making a decision about the comparative reliability of 
specific experimental data, it appeared to bo sufficient at this time to 
report the sign of the deviation. One complication, for example for 
magnesium band alloys, is that small alloy additions do not stabllisse 
the a cubic phase. Since the G/a ratio of Mg is not ideal, packing of 
individual spherical cells in which the electrons exert a pressure can 
not. be the solo criterion for determining the volume. Some modification 
of the model to include cither structure sensitive properties or perhaps 
covalent bonding may be necessary to understanding the lattice dimen- 
sions in these alloys. 

One illuminating test of this simple electron gas theory is the 
group IVb metals such as titanium, airconium and hafnium, which do exist 
in both cubic av\d hexagonal modifications. In the upper portion of this 
drawing, the a^ is plotted for the cubic phase Er-Ti alloys. However, 
in the. middle portion, two separate branches arc showi for the para- 
meter of hexagonal '/ir based and Tl based alloys. The associated c^ 
values arc presented in the lower portion of ViQj 10. Those might be 
interpreted in terms of the definitely positive “deviation from Vegard's 
Law presented in Hg. 11. Phase separation to the larger cubic phase 
would be compatible with the present model. At this time it Is not 
clear why hafnium-niobium alloys remain cubic iior why the Zr-llf do not 

form a cubic phase. 

The relative deviations 100y(G^)/p^^ can be used to reduce confusion 
about units, and values are presented in Table VII for cerium based 
alloys containing 10 per cent of several solutes. The relative devia- 
tions have the correct sign in both thorium-poor solid solutions when 
compared with the recent data of Waber, Harris and Raynor ^ but are 
smaller in magnitude than the observed values (4.14 for o’ and -1.31 
for y) . The relative deviations they computed using Friodol's Elastic 
Model^'^’^^ are similar (3.14 for cy and -0.57 for y) to those listed 

above, 

■rN.<!KUT PAGE I6a HERE. FOLLOWING END OF PAGE l6 
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Table VII 


Rolauiva Dcviauiono from Vegard's Law 
Calculated for Gerlum-liaae Alloys 


Solute 

Alpha Go 

Gamma Cc 

Magnesium 

-1.47 

i-A.75 

Aluminum 

2.07 

+4.57 

Lanthanum 

-1.84 

-0.09 

Zirconium 

H-0 . 32 

-2,29 

Tlnorium 

+2.73 

-0.89 
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Tluj sif'ti oC VO'i^) ‘loor. agireo with thou toporecd by Gschnoitlnor 

Cor dilute Mj' alLoy« but does not agree well with the Cindlngs o£ 

47 

Harris and Raynor Cor cerium dissolved in &-Zv. 

48 

Davison and Smith estimated the enthalpy o£ Cormation o£ CaMg^ 
in a relatively similar manner. They indirectly allowed Cor the shift 
in the bottom of the bond E„(R„) with concentration, The volume of 
the intermetallic compound was found from the equation 





(X8) 


Further, they assumed that was 1.225 to correspond to a Laves 

‘‘phase. In this author's opinion, it would have been preferable to 
have estimated the average number of electrons per atom and the cal- 
culated p and E. values. Tiicy used Brooks' values of R for Ca and 

c 49 ° . 

Mg. Their data are presented in Table VIII. Smith has found this 
prescription works reasonably for a number of intermetallic compounds, 


AFTER COMPLETION OF PAGE l6a. GO TO PAGE 17 
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Table VIXl. Comparison of Tlivee Physical Quantities for Laves 
Phases Calculated by fhe Prohlieh-Raimos Method 
with Kxperimcnual Values 

(After Davison and Smitl/^^) 



Heat 

of Formation 

Volume 

of Formation 

Compressibility 


-AHj. 

In kCal 


-AV 

In 

2 

X in Cm /kg 


Smith 

Brooks 

Exp. 

Calc. 

Exp, 

Calc. 

Exp. 

Ca MRg 

9.8 

11.2 

3,2 

5.66 

. 5.66 

265 

3.77 

Sr 

7.9 

9.39 

1.7 

9.64 

8.42 

— 

M M M 

Ba Mg2 

2.8 

5.69 

0.5 

15.35 

7.03 

— 

— 
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One Inut'vprcuUon o£ the falXuro oC llf-Xi* to ecparote Into a cubic 
phnao can bn made in terms of the compreaaibility of the host lattices. 
As Fig. 4 shokis, the x values for llf and gr are very similar. In con- 
trast, 'fi is significantly loss compressible, The present adaptation 
of the Frbhlieh-Ilaimos model to binary alloys permits one to calculate 
the variation of x with composition. A sot of curves for Tl-5ir alloys 
is shown ?.n Fig. 12, This result*^ rrom evaluation of Eq. and not 
from a simple scaling of the cwo P-V curves for the pure metals, We 
have not ve.rified that the experimental compressibility of alloys is as 
variable as indicated by Fig. 12. The literature indicates that very 
little data has been collected. It seems that this would be a use- 
ful test which could be made of this volume dependent theory. 

To summarise, the major fOk-tures of this FR model for a metal and 
its alloys depend on the relation between the pressure and the density 
of electrons, This in turn depends on the volume in which they are 
confined. However, any details of the spatial arrangement of the 
atoms are ignored. As such, it leads to a good interatomic force 
law. In the Frbhllch-Raimcs model it is assumed that the lowest energy 
of an electron lies between E,^, for one element and E^ for the 

other. Because of the simplicity of the FR model any Brillouin zone 
effects have not been dealt virith, 

Electron Cell Model 


In preparing this manuscript, the author came upon a very similar 

50 51 52 

recent treatment of alloys by Bolsaitis and coworkers ’ ' , It, 
too, is based on Prohlich's equations^. However, the spirit of their 
study is different since they have paramaterized all terms and in- 
cluded several further thermodynamic corrections. They represent the 
cohesive energy by the equation 


1» + + I + E 

^ R^ R^ 
s s 


AK 'i* E 
p Vib 


(19) 
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C:oc?f£iuiettu A could have beei\ obtained Irom Kq. 3 and coeCficlenu F 
could have been obtained by combininR the tomaining tenns in Kq. 3 
with the third term in Kq. 9. Since they dealt only with copper, 
silver and gold, B would have been equal to (1.2 - 0.916). Their K^^ 
is given by Kq. 14 except that the r^ and r^ were not Included in 
the numerator. They also included a van de Waals attractive term 
C/R^, and a vibrational energy term. By virtue ol: having a consider- 
able number ol' paramoter.s to adjust, excellent agreement with the 
experimental dtat can be obtained. They introduced lurthor ai'bitrar- 
inefJs when dealing with alloys by integrating the Grunoison equation 
to get an entropy correction 

C *“C T 

S' \’ — o\ w'^VB dT (20) 

They furtiher limited the ionic reviulsive term From Kq. 14 to unlike 
pairs o£ atoms. 


They varied the local charne density and cell radii according to 


vSandersou' 


'i) 


stability ratios SR and in later papers, they used a 


dcL'initLon given by Coulson for the stal)ility. Bolh definitions as- 
sumed that tlie charge densities were uniform within an Electron Coll 
and could be scaled. .Specifically, they chose Che radius of the 
copper and gold atoms acc.ox'ding to the prescription 


P (r) 
cu 

(^') 

^ au 


3 


au 


cu 


SR 


cu 


( 21 ) 


SR„ 


au 


and the moan alloy 'olurae 


V CJ 


^^cu 


u 3 , 
r 't V 
cu ^au 


0 3. 

r . ) 
au 


( 22 ) 


Since neither radius is its equilibrium value, they estimate the 
volume of Mixing by replacing r^^ and by their equilibrium 
values and r “ ii\ the parenthe.sis and then multiply It by 

C\l 

(1+Vj^j). This is an ad-hoc but plausible procedure. 
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Let us turn now to a second siraplc model which depends on knowing 
the band structure o£ each component. This model is sensitive to the 
details of tl^o crystal arrangement. In general, it is available for 
only one interatomic distance and a complicated calculation is needed 
for each atomic volume as subtle changes may occur in the relation of 
one band or K(k) curve to another. For practical reasons, the following 
model may bo regarded as independent of the atomic volume. 

LlNIiAU COMBINATION OF Dl'NSITY OF STATES 

In the basic theory of band structures, one allows one state and 

tvro spins for each atom in the piece of metal, and in general, the dis- 

23 

Crete spectrum of E(k) is replaced by a continuum if 10 atoms arc in- 
volved. One observes that when oite discusses impurity effects in semi- 
conductors, the interesting case occurs when the impurity level lies ill 

♦ 

the gap between the conduction and tlic valence bands. Since the indivi- 
dual impurities are far apart in general, and arc screened by a dielect- 
ric medium between them, the levels can be treated as discrete. However, 
if several types of defects are present, several very narrow evergy 
ranges £nd henee spikes appear in overall density of states. The height 
of these reflect to a reasonable extent how many impurity atoms of a 
given typo are prci^cnt. As a second observation, where band structure 
calculations for intcrmctallic compounds have been made, one can fre- 
quently assign energy ranges whore the states in the E(k)curvos can be 
Identified as coming primarily from states derived from one constituent. 

Examples of this are the important studies on p-CuEn alloys and isomor- 

54 55 

pliic phases by Keith Johnson and his colleagues * . The third observa^- 

56 

tioi\ comes from Sovon's study of the Coherent Potential Model ' for an 
AB alloy. This work shows a modest increase in K(E) values near ends of 
the allowed range of energies for constituent A and constituent B where 
N^CE) and Njj(E) do not overlap— o thenar i so the N(E) curves appear additive 

18 

These observations load Koskiraaki and Waber to propose as a work- 
ing hypothesis that a simple linear combination of the two density of 
states would suffice to substantially model the case of randomly distri- 
buted solid solutions. 

In order to do this adequately it is necessary to abandon the 
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tradiuional idoa of measuring all states from the bottom of the energy 
range (usually a point in k-spaee called F) . Both must be referred to a 
common energy. When the energies arc measured with respect to vacuum 

the 13,, Values for different metals, they occur in a relatively small 

fc X 

range. This Important observation ha.s becui dl.s cursed by Snow and Wf'ber 

wlieivas i.lie boLLoin of each band falls trclat l.vely) smoothly with the atom- 
ic number. The curve.s for the two pure metals are suown in Fig. 13. 


The LCDS model gives answers readily if the N(13) data is available 
for both constituents, and these alloy N(E) curves are in good agree- 
ment with tl\o experimental facts. In Fig, 14, a typical calculation of 
the composition dependent N(E) curves are presented for a solid solution 
of copper and nickel. It is very interesting that the set of these 
curves obtained with the LCDS model show many of the features exhibited 

in alloy results of Eastman'"^ ' and' Spicer”^ fsho use’d photoem'ission 

59 

spectroscopy and Clift, Curry and Thompson who used soft X-ray , 
spectroscopy. At least the morphology of the alloy N(E) is compatible 
with experience. Mo detailed comparison is ' appropriate at this stage 
of develop -.vUt. 


One observation might be made in passing. Since the individual 
values of Ej^^. are similar in energy, no large transfer of electrons to 
one typo of atom occurs. Charging effects such'as those discussed at 
the recent conference^^ will bo less than if one combined, the N(E) 
curves after the bottoms of the two bands had been aligned at the 
arbitrary' sero. 


Now to locate the Ferrni level for the alloy, one knows that the 
number of electrons is G^n^ + Cj^n^^ if there are no major charge transfer 
effects. The bottom of the alloy band w,ill occur at the lower energy of 
the two N(E) curves although it may not be heavily weighted. One finds 
the mean number of valence electrons n per atom from 
one fills up the weighted states of the linearly combined Ny^(E) and 
Ng(E) data to contain n electrons. Thus, live Ep values will have a 
definite composition dependence, as will N(Ep). The curve obtained in 
this manner for copper-nickel alloys is presented in Fig. 14. 
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For background information, tha LCDS curves of N(Ep) arc. present- 
ed in Fig. 15 for Ni and Al, and in Fig. 16 for Cu-Ni alloys. The 
decrease in N(Ep) near 60 per cent solute has been attributed to ■ 
filling of the d-states in nickel. Analysis of the data shows that 
some d-like states are unoccupied even when the concentration of solute 
exceeds the traditional amount needed at room temperature to fill 0.6 
hole in nickel. 

Examples of these N(Ep) curves are presented for several binary 
alloys below. The curve for Ti plus Al in Fig. 17, is rather uninter- 
sstitT,g since no sharp peaks and deep valleys occur in either of the 
component N(E) curves. 

The case of Ti plus V, which is shown in Fig. 18 is more inter- 
esting, particularly at the vanadium rich end. Two curves have been 
shown in Fig, 19 for the solution of Ti with Molybdenum, They are 
based on two different configurations of titanium, used as input to 
the APW program. It is more likely that the effective configuration 
in the solid BCC phase is d'^s'^ as Snow and Wabet* discuss. Despite 
whichever curve one choses to use, the shape of N(Ep) curve is surpris- 
ing since one might at first, have anticipated simple curves with 
slowly changing slopes, as one saw for aluminum additions. The reason 
for the rapid variation of • N(Ep) with composition is because the 
individual N(E) curves have a considerable amoxmt of structure. 

One even more striking curve is that presented of N(Ep) for zirconium 
and molybdenum in Fig. 20. 

Before closing this topic, another illustration of the power of 

the LCDS model might be made. Calculations were done for the BCC 

and FCC forms of several pure metals. In the case of iron, calcula- 

6 X 

tions by Koskimaki and Waber showed, that the o' or BCC form of iron 
is more stable at room temperature than austenite. This was establish- 
ed by integrating E>N(E) over the occupied states. The transition 
energy occurs with the right sign, but was slightly too large° . 
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Wa invcfiuigatci'^^ two Uypicnl alloy additiona to iron, using tlv^ 
lews model— namely, chromium and mer'ganeoe. The total alloy cnergloc 
are ohown for the two eases. The first, namely Fig, 21, chows that 
M, a well kncwit austenite former, will stabiliKC the y phase at O'^K. 
The crossover point is surprisingly close to the value determined 
from the Fo-Mn phase diagrom. In contrast, as Fig. 22 shows, no 
concentxatioii of chromium will stahili^e the gamma phase at 0*K. 


A logical extension of this work would be to study the effect of 
temperature on these two kinds of alloys. One defines formally the en- 
ergy Ep(T) the Fermi level for any temperature, as the energy at 
which the probability of a given energy state, being occupied is one- 
half-- this is more general than the usual dcfinir.ion which only per- 
tains to O'^K. Specifically, one would multiply the alloy N(E) curve 
by the Fermi distribution function 


f(E,T,Ep) = 


1 + exp 


\ KT / 


-1 


(23) 


to obtain the rcpopulation of the oloclconic levels in an alloy, With 
good luck, something approaching a phase diagram would be obtained, 


It is worthwhile to emphasize here that no thermodynamic facts 
have been built into the investigation. The atomic number and the in- 
teratomic distance are the oixly pieces of information fed into this to 
obtain individual density of states curves and the raw N(E) curves in 
Fig. 23 were then combined to get the results just presented. Inciden- 
tally,' the metals were "mixed" in the paramagnetic form and so the &- 
B Curie transformation was n'ot used lo bring about the change in 
stability of. the a and v phases^^. Many ro.iinemcnts such as a better 
treatment of exchange, and the separate handling -of spin up and spin 
down states, could bo introduced. Those would probably Lead to an im- 
proved agreement with experimental results, it is, however, striking 
that the present simple treatments outlined above do lead to useful 
predictions and increased understanding. 
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The Ik'Uinau Foynuum theorem ahowa that the Force F, acting 
on the p~ nucleua o£ a aolicl or molecule can be deduced (almost) 
classically From the total electronic charge density even though 
the latter is obtaitied From a quantum mechanical calculation. 

More speciFically, the x component oF F can be written as 


F 

px 






\’ 

q 



(R -R ) 
^ p q' 



( 24 ) 


where is the charge on that nucleus, where is the x compon 
cut of r and p(r) is the total one electron density such as 
is Found in a se3.f “Consistent calculation. That is, the torce 


on tlic nucleus is the electrostatic Force on a nucleus due to 


other nuclei and to the detailed distributimi^ oF surrounding elec- 
trons. 


While this Formulation is the adiabatic one, there have 
been dynamic treatments such as 



where W^^ is the iaternuclear conlombic Force and bV/fex^ i.s the 
instantaneous Force acting on the nucleus due to varying charge 
density p(r,t). 
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Uob' 5 ^'^ liafi wade a very up-lo-daue review of uhc force' conecpC 

* C f 

deduced from the llallwau Fcyuman iheory. Slater has rederlved Che 
Cheoroiu based on mulclplc scaecerlnp; Xo- method he has used for de- 
riving molecular wave functions. The promise of this Ilcllman-Feynman 
theory has been less than crpected and the difficulty has been attri- 
buted to the poor quality of molecular wave functions which have been 
available. 

This theory looks particularly useful for obtaining the relative 
movement and rearrangement of atoms at a free or an internal interface. 
It also appears to have a ready application to absorption of gaseous 
species on a metal. 

vSlator^^*^ pointed out that the equilibrium spacing of the surface 
nuclei will increase if there is an excess charge density on the side 
away from bulk; it will pull the nuclei out. In the contrary situa- 
tion, it will decrease the interatomic distance. The vector sum of 
foi‘ce‘3 arising from the non-spherical electron distribution v^ill give 
rise to a surface pressure. 

65 

In 1969, Wannier e_t ai calculated the foi*ce for lithium atoms 
at the surface. They found a large positive internal pressure which 
x^ould have led to considerable dilation. After pointing out some 
errors in Wannier' s treatment, Kleinman^^ showed that the forces on 
individual atoms will vanish if one assumes a slight increase of 
surface charge by 0.049e and with a slight outward movement of the 
nuclei in the surface plane by only 0.00395.i^. Without this dilation, 
Kleinman found a very substantial negative internal pressure redistri- 
bution near the surface. 
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i)LHi:USHION 


Ut 5 iu 8 ai^ ionic model, Gilwan^^ waa able Uo show that the elastic 
compressibility modulus H oi! a series o£ AB compounds was proportional 
to the inverse fourth power of r^j^. When applied to metals, it was 
found that K depended on r~'^ where v approached five. 




The compressibility is properly obtained in the following manner 


B “ 


1 

K 


C3 





12itH. 


d^'K,^, 

dtr 


(26) 


As an iilt<*rnativo procedure, one miglit consider calculating ^ from 
Equation 12 "which would load to a term in dV/dR^. However, this 
equation is already another form of 



(27) 


In either form, one would obtain an expression involving before 

the second derivative, but the numerical factor would be significantly 
differevit. 


Oiie can readily compute the second derivative of with respect 
to R„. This is 


__ 4 [_3NR^ H’ 6£(&+1) I ^ 3 1 6f (i-l-l)-4.42 a1 
dR^ R^ 


!‘3N - 1.2N -0.016 al 

H* 5 

R 

s 


4b 

f’B.M 


S. . 
. -id. 

^BM 


( 28 ) 


where a involves the cube root of N. 
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Thufi, B can be writeon as 

,2 


•12NR 


B 


12rrir 


b 2hJL(M) 


12rrU , 


^ - 13.26a ^ 

12 ., 


( 29 ) 


since both R and R are greater than unity, the lead terms become 
so 



. 3£(M'l) 

2nR^ 


<30 


ThJ.s deduction can be checked by using experimental values of B. 
Tlie slope in Fig. 24 is for various transition metals and metals v;ith 
differiiig valences and crystal structures is larger than 6. Subsequent 
plots given by Gilman^^ show that the slope is larger than 5 for trans- 
ition metals. This is merely a further indication of the type of 
svicccss which the volume -dependent FR mode.1 can be expected to yield. 


One purpose of this paper is to emphasise how directly many of the 
properties of metals and alloys are tied to the pressure of the elect)*' 
ronic gas which comes from the electrostatic repulsion between the 
electrons. Equilibrium distances, heats of fomation as well as com= 
prossibility can readily be computed using slight modifications of 
the Frbhlich-Raimes model. 

In the FR model, the kinetic energy term is computed as though 
the electrons were represented by a parabolic band, i.e. by a structure 
less N(E) curve. But we know that such a description is inadequate 
per se . To offset this deficiency, we have incorporated a discussion 
of the use of the richly structured N(E) curves for actual metals. 


Lot us turn to the LCDS model which is structure sensitive. The 

volume dependent effects are not included here, but they could readily 

be incorporated when sufficient data on the variation of the EQc) 

curves as a function of R become available. 

s 

The Rigid Band model assumes that the Fermi level and hence the 
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N(E+p.ll) anti values at Iv. arc readily calculaucd from (e/a)rat:io 

wiehouU knowing partial spin bands o£ the constituents, 

One successCuI aspect of the LCDS siodel is that it offers a 
simple logical reason for the deviations from the idealited Slator- 
Pauling curve of unpaired spins N(E). For differing atomic numbers, 
the N(E) curve may bo very different in width and variable in location 
with respect to vacuum. The number of paired spins at Ep is thus 
not simply related to the density of states curve of either constituent 
since It depends on the response of the two partial spin up and two 
spin do\m bands in a binary alloy to the internal magnetic field pH. 
Even in a pure metal, N(F,pt) is not equal to NCEpI) and hence 
transfer of electrons from one partial bemd to another occurs and 
leads to an unbalance in the pairing of spit\s. If the two constituents 
arc elements of very different atomic ttumber, the effect of a magnetic- 
field on each N(E) curve may be quite different while the affect of 
solute additions may result in recogiaiaable trends in Slater-Pauling 
curves. There is no reason to anticipate that the height and the width 
of the ii\dividual alloy curves could be superimposed when plotted 
agaiirst (e/a) values. 

R. E. Hummel et , usiiig differential reflectivity, obser-ved 
that Che adsorption edge (or change in the imaginary component c., of 
the dielectric constant) associated with the d-bands of copper, did 
not move to lower energy with nickel additions. The energy separation 
octwoeu the peak energy did ivot move, Tliey also measured the reflec- 
tivity of GU"Zn alloys and noted that the energy soparatioii between 
the d-baiids and Ep does increase with Zn content. This they rational- 
itc by making the energy of the state the same in pure C,u as in 
o'-brass. Thus, the energy level Ep should increase with due to 
the increase in the mean number of valence electrons. These observa- 
tions on CuNi and GuZn alloys arc not inconsistent with the LCDS model. 

Other theoretical predictions for CuNi alloys have been made by 
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Lang and Ehronrcich^^ and by Stocks, Williams and I’aulkner^^, The 
former group predicts that N(Ep) should decrease linearly with 
whereas the Oak Ridge group which used the Gl»A approach predict a more 
rapid initial decrease followed by a leveling off above 60 per cent 
copper. It is interesting to note that lluffner cj; ad' ^ point out that 
their experimental density of states for the alloys "...to a very good 
approximation, be made up by aupcriimposing those of Ni and Cu." 

To date, no synthesis of these distinct approaches has been 
attempted. Pragmatically, they are two limiting cases of the more 
complex true situation. While my colleagues and I are working on such 
a synthesis, it is hoped that this review will stimulate some readers 
to understake the overall problem, and make a significant contribution 
to the science of metals. 
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Fl8. 1. 


Fig. 2. 
Fig. 3. 
^ Fig. 4. 
Fig. 5. 


Fig. 6. 


Fig. 7. 


Fig. 8. 
Fig. 9, 


Plot illustrating Che dependence ei U on the ioniaation 
2 

potential w and on the valence N. The continuous curve 
is band N«4, Waber and Larson used r^ in place of 
Their Fq. (7) corresponds to Eq. (6) herein and their 
Eq. (5) is obtained by setting A^O. 

Plot of the different electron-electron interaction toms 

7 

in Eq. (9) (after Larson and Waber ). 

Determination of and p from E^ and E.p curves of 
tiianium and airconlum. 

The calculated pressure dependence of x fo^^ four different 
tetravalont metals. 

Comparison of experimental shock wave data of McOueen and 
Marsl/*^^ Cor nickel with two calculated P versus V curves. 
The higher valence N**3 (wliich is chemically stable) is on 
approximate way of taking account of repulsion between 
d-olectrons and ion-cores. 

Similar Shock-Wave data for gold compared with those 

calculated from K„,(R ) with and without E.^„. from 
T 8 ionic 

Eq, (]4) included. 

% 

Plot of evaluated at the composition dependent 

radius P for titanium-cerium alloys. The effect of 
ns.suming that cerium is tetravalenU is to cause only a 
slight depression in the E.j, curve. It is interesting that 
p (C^) lies above the straight line connecting the p 
values for the two pure materials. A positive deviation 
is indicated from Vegard's Law of additive atomic radii. 
Similar plot for sodium-silver alloys which indicates a 
large' positive deviation from Vegard's Law. 

The .strong negative deviation from Vegard's Law for sodium- 
magnesium alloys contrasts with the alloys curves presented 
so far, 
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Pis. 10 


Fig. U. 
Fig. 12. 

Fig. 13, 


Fig. 1^. 


Fig. 1';. 


Fig. 16. 


Fig, 17. 


Pig. 18. 
Fig. 19. 


Fig. 20. 
Fig. 2'i. 


M 


E.sporimental valuct; of chc two Xaeticc pararactcrts of a 
series of hcxagooal mucalo, In thu middle of tlic com- 
poettlon vange, titaniutn-zircoitium alloyo become body 
eenteirod cubic au room temporature according to Duvea 
Plot of and p <C^) for titanium-zirconium alloys. 

Variation of the pressure volume curves calculated for 
ti u an ium- zirconium alloys . 

N(K) curves for pure copper and nickel obtained by Waber 
and Snow^^ using Slater's exchange parameter xff»l. Those 
two curves arc combined to obtain Fig. 14, 

Set of density of states NCE^,), curve.s calcul.ated for several 
nickel-copper alloys using the LCDS model of Koskimaki and 

,, 1 la 

Waber . 

Composition dependent change in the density of states 
at the calculated Fermi level for homogeneous 
nickel -aluminum alloys, 

A similar plot of the variation of N(Ep) for homogeneous 
nickel-copper alloys. The individual Fermi levels are 
not shown in Fig, 14. 

Composition dependent density of states at the Fermi 
level Ep calculated for homogeneous titanium-aluminum 
.olloys. 


.Similar plot of N*Ep)for titanium-vanadium alloys. 

Similar plot of N(Kf,') for titanium-molybdenum alloys. 

*■' 2 2 
Separate curves are given for the configurations d s and 

o ]4 

d s assumed for BCC titanium by Waber and Snow' , 

Similar plot for homogeiieou.*? zirconium-molybdenum alloys. 
Total energy of occupied electron states at O'k as a 
function of composition using bonds calculated for both 


, 6 / _ 2 , 


2. 


pure phases with iron (3d 4s ) and manganese (3d“4s ). 
The transformation from FCC to BCG as manganese is added 
is Indicated. 
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Pig. 22, A iiimilar plou which shov;y that, no chror,»iuin addition ic 

ouiWclenU to ctablUzo the auotmiitc at 0 "k . Based on 

« # 

data calculated by Snow and Kabet^ ^ using chmniutn 
with the ? input configutaulon. 

Pig. 23. Compariso.- of the density of GUates curves for FCC and 
BCC iron in a paramagnetic state. Note the difference 
in height and slope at Ep. (after Snow ond Wobor ) 

Fig. 24. Log-log plot of experimental values of the bulk modulus 
ag'ainct the cell radius r^ (in Angstroms). Slope of 
the line is larger than 6, 
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